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A recently proposed scheme is used to saturate the spectral side of the QCD sum rules derived 
from the thermal, two-point correlation functions of the vector and the axial-vector currents. At low 
temperature, it constructs the spectral representation from all the one-loop Feynman diagrams for 
the two-point functions. The old saturation scheme treats incorrectly some of these contributions. 
We end up with the familiar QCD sum rules obtained from the difference of the two corresponding 
vacuum correlation functions. The possibility of obtaining new sum rules in other media is discussed. 

I. INTRODUCTION 

In their original work extending the vacuum QCD sum rules M to those at finite temperature, Bochkarev and 
Shaposhnikov recognised the importance of multiparticle (branch cuts) contributions in addition to those of single 
particles (poles) in constructing the spectral representation of thermal two-point functions. Thus for the vector current 
correlation function, they included not only the p-pole with temperature dependent residue and position, but also the 
7T7T continuum. In the same way the spectral representation for the nucleon current correlation function would consist 
of the nucleon pole with temperature dependent parameters and the 7riV continuum. Although such a saturation 
scheme is quite suggestive and had been extensively used in the past HQ], it lacks a theoretical basis, leaving one to 
suspect if some equally important contributions are left out. 

A definitive saturation scheme emerged with the work of Leutwyler and Smilga ||, who calculated the nucleon 
current correlation function in chiral perturbation theory. Being interested in the shifts of the nucleon pole parameters 
at low temperature, they considered all the one-loop Feynman diagrams for the correlation function and evaluated 
them in the vicinity of the nucleon pole. Koike Q examined the contributions of these diagrams in the context of 
QCD sum rules. He found a new contribution arising from the nucleon self-energy diagram to the sum rules, not 
required in the saturation scheme mentioned above. 

The one- loop Feynman diagrams for the nucleon correlation function were further analysed in Rcf. 0. In this set 
of diagrams for the correlation function, one has not only diagrams with the (single particle) pole alone and the (two 
particle) branch point alone, but also other (one particle reducible) diagrams, appearing as product of factors with the 
pole and the branch cut. As an example, take the case of a vertex correction diagram having this product structure. 
It may be expressed as the sum of the pole term with constant residue and a remainder, regular at the pole. Clearly to 
find the pole parameters one may confine oneself to the pole term alone. But if one wishes to evaluate the diagram for 
large space-like momenta - the region of relevance for the QCD sum rules - both the pole and the remainder become 
of comparable magnitude. It is these remainder terms which are not included in the earlier saturation scheme. 

In this paper we write down the thermal QCD sum rules following from the vector current and the axial vector 
current correlation functions, constructing the spectral side from the set of all one-loop Feynman diagrams. At low 
temperature only the distribution function of the pions is significant in the heat bath. Thus, of the two particles in 
the loop, at least one must be a pion, the other being any one of the strongly interacting particles with appropriate 
quantum numbers in the low mass region. The vertices occurring in the diagrams are obtainable from the chiral 
perturbation theory ||,|9| as well as from other formulations of the effective theory of QCD at low energy |[o]]. But 
the former theory alone is based only on the chiral symmetry of QCD, the others assuming one or other unproven 
symmetries or relations. Thus the chiral perturbation theory is singled out to determine the vertices, if we claim the 
sum rules to follow from QCD. 

As with the nucleon sum rules Q, we subtract out the vacuum sum rules from the corresponding full sum rules 
at finite temperature, equating, in effect, terms of 0(T 2 ) on both sides. All our calculations are done in the chiral 
symmetry limit, though we keep nonvanishing pion mass in intermediate steps. 

In Sec. II we collect some results to be used later. In Sec. Ill we analyse the one-loop Feynman diagrams to 
construct the spectral representation for the correlation functions. In Sec. IV we use the known results of Operator 
Product Expansion and write the sum rules. Finally our concluding remarks are contained in Sec.V. 
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II. PRELIMINARIES 



Here we review briefly the kinematics of the two point correlation functions of the vector and the axial-vector 
currents. Then we write the interaction vertices from the chiral perturbation theory and the results of calculation of 
a prototype loop integral. 



A. Kinematics 



Consider the two thermal correlation functions, 



T^ = iJ^xe i ^TrgTV^(x)V^O) , (2.1) 

and 

T^ b = i J d 4 xe^TrgTA«(x)A b u (0) , (2.2) 
of the vector and axial vector currents, 

T a a 

Y?(x) = q(x) ltl —q(x), A*{x) = q(xh^ 5 — q(x) , 

generated by the SU(2) flavour symmetry group of the QCD Lagrangian. Here r a are the Pauli matrices and 
q = er^ H /Tr e~^ H is the thermal density matrix of QCD at temperature T — 1/(5. Note that in the limit of chiral 
symmetry, in which we shall work, the axial vector current is also conserved and the kinematics, in particular, the 
invariant decomposition is the same for both the correlation functions. 
The current conservation leads to the invariant decomposition 

T^M = S ab (P^T t + Q^Tt), (2.3) 

where the gauge invariant tensors are chosen as 

p .<h^_<f_~ . n _<f_~ . 

qz qz qz 

with Up = — ujq^/q 2 , where u M is the four-velocity of the medium and u and q are Lorentz invariant scalars, 
u> = u ■ q and q = \Ju) 2 — q 2 , representing the time and space components of q^ in the rest frame of the heat 
bath (uq = 1, u = 0). The invariant amplitudes are functions of the scalar variables, say, q 2 and u>. They can be 
conveniently extracted from the Feynman diagrams by forming the scalars, 

Ti=g» v T^ , T 2 =u' i u u T llv , (2.4) 

which are simply related to the invariant amplitudes. 

Having cast the kinematics in a Lorentz invariant form, we choose to do calculations in the rest frame of the heat 
bath. The kinematic decomposition (2.3) leads to a constraint on the invariant amplitudes which, in this frame reads 
as 

T t (q , q = 0) = q 2 T^q = 0) . (2.5) 
Using this equation, the two sets of amplitudes, T/ jt and T^a can be related for q = as 

Ti = ^T 2 , T t = -iTi . (2.6) 
q z 3 

Note also the symmetry of the imaginary parts of the amplitudes, 

ImT lit (-q , q = 0) = ImT lit (q , q*= 0). (2.7) 

In the real time thermal field theory, we are going to use here, each of the above amplitudes stands for a 2 x 2 matrix, 
whose components depend on a single analytic function |]TTf . This function, in turn, is determined completely by the 
11-component function itself: their real parts are equal and the imaginary part of the former equals 7r tanh(j3qo /2) 
times that of the latter. (The factor 7r _1 is inserted for convenience.) To avoid further symbols, we henceforth redefine 
T to denote this analytic function. Its spectral representation at fixed q is given by, 

T t A^= r dq '° 2 '' mT ^ ,qi - 

Jo % 
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B. Dynamics 



As we pointed out already, we must describe interaction among particles according to the chiral perturbation theory. 
The calculation of the two point functions in this theory is most conveniently carried out in the external field method 
PL where the original QCD Lagrangian is extended by introducing external fields v^(x) and a°(i) coupled to the 
currents V°(x) and A°(x), 

Cqcd - C QCD + v a Jx)V»{x) + a a Jx)AZ(x). 



The resulting interaction vertices are obtained in Ref. [g_2||13j for the symmetry group SU(3)r x SU(3)l- 

At low temperature the pions dominate the heat bath. We thus consider the reduced symmetry SU{2)r x SU(2)l 
and write down the chiral couplings of pions with the external fields and the observed particles, to be encountered 
in the one-loop diagrams fl4j| . We leave out the interaction among the pions themselves, as it would give rise to 
corrections of order higher than T 2 . Their interaction with external fields are given by 

£; nt (7r)=A,(7r)+£ (7r) (2.9) 

where 

C v (n) = ■ TV X d^TT + — (tV ■ TV V ^ ■ — V ^ ■ TV ■ Tv) 

C a {-K) = -F^ap ■ a* - F^a^ ■ d^Tv - ~ (tv ■ tv ■ - a M • tv a 1 " ■ tv) 



— {TVTvd^TV-a^-TV-d^TVTVa^ , (2.10) 



the letters in bold face denoting isospin vectors. The third and the fourth terms in C a (Tr) represent corrections 
respectively to the first and the second term, when two of the pion fields without derivatives are contracted in them. 
Such a contraction can be worked out from the thermal pion propagator 

A n(fc) = T^~^ + 2 ^( fc2 " ml)n(k) : (2.11) 
h — m% 

where n(fco) — (e^ ' — 1) _1 is the Bose distribution function. Then the thermal part of the contraction gives 

TrQTTV a (x)TV b (x)\u^S ab J -0^n(k o )S(k 2 - ml) = 5 ab ^, (2.12) 

in the chiral limit. We thus obtain the effective form of C a (iv) as 

4.M = ^ (l - ^) »n >a"-F v (l- a, ■ &>* . (2.13) 

Next we write the couplings of the isotriplets [p(770), ai(1230)] and isosinglets [w(782), /i(1282)] of vector (1~~) and 
axial vector (1 ++ ) mesons respectively. Although we take their fields to transform according to SU(2) in constructing 
their interactions, we take the physical (zero temperature) masses of the multiplets of each of the SU (3) octets to be 
degenerate. Then the coupling linear in the vector meson fields are given by 

^2 



1 - Wv" ■ [d M p v - d vPv ) + j-dW ■ [d» Pv - d vP)1 ) x tt| 

■^rf^ ■ x U tv - y^Le^^Tv ■ d*v°, (2.14) 



while those linear in the axial vector meson fields are 

C(A) = -H*L ((l- d»a v ■ (d^a lu - d v a llM ) + -L&>v» ■ (d^a lu - d u a lfi ) x tv 



■mAr-n 



e^ a tfd»n-d x a°. (2.15) 



3 



Finally the quadratic couplings of the triplets with the singlets and between themselves are given by 



C(V, A) = -2e^ ACT ( §-d»u> u p x ■ 5 ct tt + f-^/jV ■ d a iz ) + f-c>V 'Kx d„7r - oi v x d^). (2.16) 



The coupling constants in the above interaction terms can be determined from the decay rates of the particles f^,[14| . 
The propagator for the heavy meson can be taken to be the one in vacuum. 

There also appear two other spin-one mesons /ii(1170) and &i(1235) having quantum numbers 1 H Jig] , But P and 
C invariance forbids their appearance in our diagrams E§. The pseudoscalar SU{2) singlets r/(547) and r/(985) also 
have no couplings relevant for us. The scalars ct(400 — 1200) and ao(980) do have vertices for our diagrams, but they 
do not contribute to 0(T 2 ), as do many of the vertices written above. We shall omit the vertices with the scalars 
altogether from our discussions to follow. 

With the above interaction vertices we can write the Feynman amplitudes for the correlation functions. A loop in 
the diagrams makes a contribution of the form 

F^(Q)=ic J -0^UAq,k)A[l\k)A^\q~ k), (2.17) 

where the particle A can be a heavy (spin one) meson or the pion itself. The tensor is given by the interaction 

vertices and any tensor structure in A u . Out of the interaction vertices we have isolated the coupling constants in 
c. Being gauge invariant, Fp V allows one to construct the invariant amplitudes Fi.t in the same way as we did above 
forT^. 

As we shall see below, the one-particle irreducible diagrams consist only of loop integrals like (2.17), while the 
reducible ones are given by such integrals multiplied with a heavy particle pole (of first or second order). In the 
former case we may evaluate directly the leading contribution of the integral and then Borel transform. In the latter 
case it is more systematic to use the exact spectral representation of the loop integral. We then Borel transform the 
resulting complete amplitude and extract its leading term. 

Until now we have considered full amplitudes, including the corresponding vacuum amplitudes. Our thermal sum 
rules need only their temperature dependent parts. So we evaluate only these (convergent) parts of the loop integrals, 
but continue to denote them by the same symbols as for the full amplitudes. 

Consider first the case where A is a heavy particle of mass ran ■ Let the full amplitude be given by just the loop 
integral (2.17). Going over to the invariant amplitudes, its thermal part is as 

r, , s f dik S(k 2 - mi)n(k) „ . , 

F ^ - ~ c j (2^ (g -fc) a -4 /M(g ' fc) ' (2 - 18) 



whose leading behaviour for large space-like momenta with q = (q 2 . — E 2 = — Q 2 < 0) is given by {ujk = y k 2 + m 2 ) 



2c f d 3 kn(k) 
+ m 2 H J (27r) 3 2uV 



^ t (0 2 ,? = o)- 7 ^— -r I TK^y-fiAQAkl)- (2-19) 



It is clear that only if fi,t(Q, \k\) is constant in k, is the integral of order T 2 ; otherwise, it is of higher order. 

To obtain the spectral representation of the loop integral for q — 0, we note that the cuts in E 2 plane in this case 
are given by < E 2 < (mjj — m^) 2 and E 2 > (ran + m^) 2 ■ The imaginary parts across both these cuts are given by 



\/ LO 2 — TTL 2 

ImFi^E) = 87T 2 E w ) ( 2 - 2 °) 

where uj = (E 2 — m 2 H + m 2 )/2E. We now have the spectral representation for F/. f given by an equation like (2.8), 
the integration running over the two cuts stated above. 

If the particle A is also a pion, the leading term of the loop integral may again be found in the same way as we did 
above for the heavy particle. However, to find the spectral representation for q — 0, one must begin with q ^ and 
then go to the limit. This evaluation is done in the Appendix. 
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III. SPECTRAL REPRESENTATION 



Here we analyse the different one-loop Feynman diagrams for the correlation functions and extract the leading 
thermal contributions (to order T 2 ) to the spectral side of the sum rules. These diagrams are the same as those 
considered in JlJ] to find the shifts in the pole parameters of the p and a\ mesons. The difference lies in the evaluation 
of the diagrams: While we evaluated them earlier in the neighbourhood of the respective poles, now we have to do so 
at large space-like momenta. 

The diagrams for the correlation functions can be grouped into three types, namely those with intermediate states, 
vertex corrections and self-energies. To write the Feynman amplitudes we anticipate the following gauge invariant 
tensors, 

A^q) = -g^u + q^qu/q 2 = + Q^/q 2 , 

B^(q, k) = q 2 k ll k v - q ■ fc(g r Al fe v + k^) + (q ■ kfg^, 

C^(q, k) = q A k ll k v - q 2 (q ■ k^q^K + k^q v ) + (q ■ k) 2 q^q v . (3.1) 

We now consider the diagrams separately for the correlation functions of the vector and the axial vector currents. 
Although we need only the T-dependent parts of the diagrams to order T 2 , we shall write the pole amplitudes in full. 



A. Vector current 

Here the pole and the one-loop diagrams are shown in Figs. 1-4. 



t\r\r\r\i u\s\r\r\ ryru^AU - i/xrvnn o/vataj- 





— 



(a) (b) ( C ) 

FIG. 1. p pole and constant vertex correction 



a, a,, % 

(a) (b) 
FIG. 2. Intermediate state diagrams 



p 

to a,, it wvrvr^ '\r\r\r\j — m a , „ 
(a) (b) 

FIG. 3. Vertex correction diagrams 



to, a, , n 

(a) « 

FIG. 4. Self-energy diagrams 
The p pole amplitude with its constant vertex corrections (Fig. 1) is given by 

T$ (q) = ~ ( ( 1 - ^s) (3.2) 

The amplitudes for the intermediate state diagrams (Fig. 2) are generally of the form of Eq.(2.17). In the case of 
the 7T7T intermediate state, one has to include also the tadpole diagram of Fig. 2(b) giving 

T^\q)=i J 0^{(2k-q)»(2k-q) u A$(k)A^(q-k)-2ig^A$(k)}, (3.3) 
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whose thermal part is 



r™w n , , ?, f d4k 8(k 2 - m 2 )n(k) T 2 , s 

T ^= V(-^)y (2^)3 ( g - fc) 2 - to 2 ~" 6~ ^ ^ 

at large space-like momenta. For the nai intermediate state, the amplitude is given by Eq. (2.17) with 

c = 2(F ai /m A F^) 2 , /„„ = q 2 (q 2 - 2q ■ k)A^ - B^, 

getting 

fl,t = E 2 (l, E 2 ) . 

Then Eq. (2.19) gives 

\m A ) Q 2 + m 2 A ' 6F* 1 ' 

The (Borel transforms of the) amplitudes given by Eq. (3.2,3.4,3.5) constitute all of the contributions to the spectral 
side of the sum rules. In the rest of this subsection we verify that none of the remaining diagrams with iru) intermediate 
state, vertex corrections and self-energies contribute to order T 2 to the sum rules. Clearly it suffices to show this 
behaviour for any one of the invariant amplitudes, say Ti and we omit its subscript in the following. 
For the ttlo intermediate state, the amplitude is again given by Eq. (2.17) with 

c^^iH^/mvF^) 2 , f^ = q 2 k 2 A^ + B^ v , f = -2|fc| 2 /3 , (3.6) 

where and below / stands for /j. Then Eq. (2.21) shows immediately that this amplitude is of order T 4 . 
Considering next the vertex corrections of Fig. 3, each of the amplitudes is of the form 

F q 2 

9 t - ^^^ , (3-7) 
uiy q z — TUy 

where A M „ is again a one-loop integral of the form of Eq. (2.17). Consider now 

E 2 



E 2 - 



m 



ME 2 ), (3.8) 



where the invariant amplitude A(E 2 ) satisfies the spectral representation 

2 f ImA(E')dE' 2 

HE)- J E , 2 _ E2 , 

which may be used to write the expression (3.8) as 

m 2 / A(m 2 / ) f dE' 2 E' 2 Imk(E') 



E 2 -m 2 v J (E' 2 - my)(E' 2 — E 2 ) 



(3.9) 



separating the pole term from the one regular in the vicinity of the pole. On the other hand, if we go to large space-like 
E 2 = —Q 2 and take the Borel transform, it becomes 

e -m 2 v /M 2 rdE 2 Imk{E) 



M 2 



I V - K J) (-™ V + E>e-^-^) , (3.10) 



where M 2 is the Borel variable replacing Q 2 . 

Eqs. (3.9) and (3.10) allow us to compare how the same amplitude behaves near the pole and at large space- 
like momenta in the form of the Borel transform. Consider the vertex correction from the ttw loop for which c = 
— (i^H^gi/mvF 2 ) and the expressions for f^ u and / are identical to those in Eqs. (3.6). Then we see that in 
Eq. (3.9) A(m v ) ~ 0(T 4 ) and the second term is indeed finite at E 2 = m v (and of order T 2 ). Next consider the 
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behaviour of the Borel transform (3.10). Because the branch point of the function A(E 2 ) also starts at E 2 — m v , its 
leading term is given by 

-my /M 2 



l-^j^^Jd&ImME), (3.11) 

which is of order T A . Thus to 0(T 2 ), the vertex correction diagram with the ituj loop contributes neither to the 
pole parameters nor to the Borel transform to 0(T 2 ). Notice, however, that while near the pole the second term in 
Eq. (3.9) can be ignored, for large space-like momenta both the terms assume equal importance. 
Consider next the nai loop with 

c = -(4F ai g 3 /m A F%), f„ v = q 2 q ■ kA„ v + B^, f = E\k\ . 

Again we see that at E 2 = m v , the residue A(m^) ~ 0(T 4 ). Here the branch point of A(i? 2 )is is at E 2 = m\. So 
the Borel transform is clearly of order T 4 . As shown in the Appendix, the irir loop contribution is also at least of 
order T 4 to the pole residue and Borel transform. 

Finally the self-energy diagrams of Fig. 4 can be analysed in a similar way. Each of the diagrams contributes an 
amplitude of the form 

F \ 2 o 4 

— ' T^-^n^te) • ( 3 - 12 ) 

my J (q z - m v Y 

The self energy is again of the form of Eq. (2.17). Using the spectral representation for its invariant amplitudes, 
we get 

C> H(E) 



(E 2 - m 2 yf 



i v f dE' 2 ImU(E') t m 2 v f dE' 2 (2E' 2 - m 2 v )ImU(E') 



(E 2 - ra 2 v ) 2 J E' 2 - m% ' E 2 - m 2 v J {E' 2 - m 2 v ) 2 

dE' 2 E' 4 Imn(E') 
(E' 2 — my) 2 (E' 2 — E 2 ) ' 



(3.13) 



which is the appropriate expression to study the neighbourhood of the meson pole. If we now go to space-like momenta 
and take the Borel transform, we get 



^ J dE 2 Imn(E) 



E 1 \ ml ml _ 1 e - mv/M 2 + ( E 2 \ e -E*/M* 



M 2 E 2 



(3.14) 



Consider first the ttlv self-energy loop, whose evaluation reveals the inadequacy of the old saturation scheme. For 
this loop we have c = 4(<7i/F 7r ) 2 and the expressions for and / are again given by Eqs. (3.6). Near the meson 
pole, the three integrals in (3.13) are all finite and are of order T 4 ,T 2 and T 2 respectively. Thus the pole position 
does not shift, but the residue does to order T 2 [Q. On the other hand the Borel transform (3.14), which simplifies 
to leading order as 

is clearly of order T 4 . To summarise, the self-energy diagram with the ttoj loop gives a correction to the pole residue 
to order T 2 , but the Borel transform of the full amplitude has no contribution to this order. By contrast, the old 
saturation scheme would retain the correction to the residue also in the Borel transform. It amounts to the neglect of 
the third, regular term in Eq. (3.13), which is of course justified near the pole, but not for large space-like momenta, 
where the pole and the regular terms are not only of comparable magnitudes, but actually cancel each other in the 
leading order. 

The nai self-energy loop has 

c-2(.g 3 /F0 2 , /„„ = - C^/m 2 Al f = -|fc| 2 (2 + E 2 /m 2 A )/3 . 

We see that this amplitude contributes neither to the pole position nor to the Borel transform to 0(T 2 ). The same 
is the case with the tttt loop as shown in the Appendix. 
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B. Axial- vector current 



Here the Feynman diagrams are shown in Figs. 5-9. The new feature here is the existence of the pion pole in 
addition to the one for the axial vector meson a\. 




FIG. 5. 7r pole and constant vertex corrections 




(a) (b) ( C ) 

FIG. 6. ai pole and constant vertex corrections 
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FIG. 7. Intermediate state diagrams 




FIG. 8. Vertex correction diagrams 



... ,-. a i f Vi .-. .-. .-, 



FIG. 9. Self-energy diagrams 
The 7r-pole with vertex corrections of Fig. 5 is given by 



(i/r ) 



2 

' (3.15) 



in the chiral limit. The ai-pole with vertex corrections of Fig. 6 is given by 



-V • ( 3 - 16 ) 



\rriAj V" 6F 2 Jq 2 -m 2 A ^ ' 
The contribution of the 7rp intermediate state of Fig. 7 to the invariant amplitudes T^ t are given by 

h * ~\m v ) 6F* Q 2 + m 2 v - [6 <} 

As in the case of the vector current correlation function, one can show that none of the remaining diagrams contribute 
to order T 2 at large space-like momenta. 



IV. SUM RULES 



We now turn to the operator side of the sum rules. Up to dimension six, the Lorentz scalar operators appearing in 
the short distance expansion of the product of two currents are 1, m q qq, G^ U G^ V and two four-quark operators. Of 
the latter two, one arises from the Taylor expansion of the quark bilocal operators arising in the Wick expansion of the 
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currents. This piece turns out to be isospin scalar. The other four-quark operator arises from the Feynman diagram 
to second order in QCD perturbation expansion, where the large external momentum flows through the internal gluon 
line. 

In the medium there are additional operators, which are Lorentz non-scalars in the rest frame of the heat bath. 
They can, of course, be written as Lorentz scalars in a general frame by using the four- velocity vector u^. They first 
appear as dimension four, which are qi/iq, u^u^Q^ and u* l u v ®9, where 0^ are the energy momentum tensors of the 
quark and gluon respectively. More such operators arise at dimension five and six, but generally contain derivatives. 

In the process of subtracting out the vacuum sum rule to extract terms of order T 2 , the unit operator drops 
out. Also drops out m q qq in the chiral limit. The isospin scalar four-quark operator cannot have a temperature 
dependence in the chiral limit. The operator qi^q does not contribute, as we do not have non-zero chemical potential 
for the fermions. The thermal expectation value of G /J1/ G A " / , 0^„ and 0^„ are all of order T 4 . Also the expectation 
values of other operators with derivative(s) are of higher order than T 2 , as each derivative gives rise to an additional 
power in pion momentum. 

Thus we are left with only the isospin non-scalar four quark operator as the relevant one for our sum rules. Their 
Wilson coefficients are known |jj , 

i J d^xe^ (TV«(x) V»(0)) - 8" b {- 9lxv + ^ ^L(O a ) 

i j 'd'xe^iTA^Aim^d^t-g^+^j^iOy) (4.1) 

where 



O v = otsq-y^ — ^-qq^ — ^-q 
A = a s q-f fJ ,j 5 — -qqj ti j 5 — Y q 



and a s = g 2 /47r is the strong interaction fine structure constant. 
The thermal average of any operator may be expanded as 

(O) = (0|O|0) + £ J l^- (n a (k)\0\n a (k)) + • • • , (4.2) 

where the sum is over the isospin indices of the pion. The pion matrix elements of Oy.A are easily worked out using 
PCAC and Current Algebra to get, 

{O v ,a) = (0|Oy, A |0) T ^(0|Oy - O A \0) 

It is now simple to write the sum rules for a correlation function by equating the Borel transform of the spectral 
representation to the operator product expansion. It turns out that both the vector and the axial vector correlation 
functions give rise to the same sum rules 



F 2 e -m* v /M> + F 2 ie -ml/M> + p 2 = _-^_(0|O y - O A \0) 



-mlFy™v/M* + m a j^e-mi/M- = ^(0\O v - O A \0) (4.3) 

It is interesting to observe that these sum rules are nothing but the vacuum sum rules derived from the difference 
of the correlation functions of the vector and the axial- vector currents, whose spectral side consists of the pole terms 
from 7r, p and a\ exchanges. Also as M 2 tends to infinity, we recover the well known sum rules 

F 2 P - Fl ~F 2 = 
m 2 v F 2 - m 2 A F 2 ai - (4.4) 

originally derived by Weinberg p7[ as superconvergence sum rules for the difference of the spectral functions. 
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V. CONCLUDING REMARKS 



In this work we have obtained the sum rules following from the two point functions of the vector current and the 
axial vector current at finite temperature. For a reliable estimate, we calculate their spectral sides using the chiral 
perturbation theory. The resulting Feynman diagrams (to one loop) can be readily evaluated for their leading thermal 
contributions. 

Though there is a multitude of diagrams to begin with, only a few of these actually contribute to order T 2 . This 
is due to the presence of derivatives on the pion fields in the interaction vertices derived from the chiral perturbation 
theory. Among the contributing ones, there is a one-particle reducible, self-energy diagram, whose contribution at 
space-like momenta is zero to order T 2 . But the old saturation scheme would derive from this diagram a shift to this 
order in the residue of the meson pole. 

The older strategy of determining the temperature dependence of the pole parameters from such sum rules is not 
relevant anymore, as we have already included the diagrams responsible for this dependence in our calculation of 
the spectral functions. Our evaluation reproduces the sum rules which follow from the vacuum correlation functions 
themselves. 

The present saturation scheme has been used also for the nucleon correlation function at finite temperature, again 
reproducing the results obtainable from the corresponding vacuum correlation function [HfjJ. Thus we now have 
enough confirmation of the correctness of the QCD sum rules in a medium. 

Can we get new results from the sum rules in a medium? It would seem that we are merely reproducing the old 
results in a more complicated way: What one could have from single particle exchange diagrams in vacuum, are now 
derived from one loop diagrams at finite temperature. But we do not believe the situation to be so. Indeed, having 
been assured of the correct procedure to saturate the spectral side, we may apply the sum rule technique to other 
media, like the nuclear medium by introducing the nucleon chemical potential. There are no reliable estimates of 
higher dimension operators like the four-quark operators in such media. It is these quantities which should be readily 
available from such sum rules. 



APPENDIX 

Some care is necessary in constructing the spectral representation for the tttt loop for q = 0. The loop integral gives 
rise to cuts in the q 2 plane for arbitrary q extending over < q 2 < \q\ 2 (short cut) and q 2 > 4(m 2 + \q\ 2 ) (unitary 
cut). As q — > the short cut shrinks to a point. But the spectral function can be singular in this limit, so that the 
part of the spectral representation due to this cut may be non- vanishing. 

Consider first the tttt loop in the vertex diagram. It gives an amplitude of the form of Eq. (3.7), where 

A^(g)=ic J ^UAq,k)dg(k)A$(q-k), (A.l) 

with 

c = (4G p F p /m 2 v F 2 ) , / M „ = (^ - 2k^)(q 2 k u - q ■ kq v ) . 

As in Sec. II the imaginary part of the amplitude may again be obtained directly by integrating over ko- But, unlike 
the case of loops with one pion and one heavy particle, the angular integration in the resulting integral over k gives 
rise to a 9- function constraining the limits of |fe|. The imaginary part of the invariant amplitudes may then be worked 
out to give 

Imki \ 1 f +1 J (~q 2 x 2 l2\(^ 2 \ 2 n fKn . 

and 

JmA, \ _J_ f°° ( -q 2 x 2 /2 \ ( 1 1 \ 2 

Imh t ) 8(2tt) 2 J l \ -q 4 /6 J y e ^\S\x-q )/2 _ x e P(\q\x+q )/2 _ i J ' q <U ' 

On the unitary cut the cc-integrals are thus finite as q — > and because of q 2 (q A ) in the integrands, their Borel 
transforms will have no contributions to 0(T 2 ). But on the short cut the ^-integrals diverge as q — > 0. Let us then 
consider the full Borel transformed amplitude with q^O, 

M - 1 ^dqle-^H^A • (A.4) 



A* ; " M 2 L " yo I ImA 
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Defining new variables A and u by qo — X\q\ and \q\x = u, it becomes 



where 



A t J 8(2tt) 2 M 2 



D (A.5) 



D= 1 f 1 



\\q\ ^ e /3(u — A|gH)/2 _ I e /3{u+\\q\)/2 _l 

_ d / l 



du V e' 3 "/ 2 - 1 7 ' 



(A.6) 



as |<f| — > 0. Thus the short cut also cannot contribute for q = 0. 

The 7T7T loop in the self-energy correction is given by Eq.(3.12) where n M „ is of the form of Eq.(A.l) with c = 
(2Gp/myF%) 2 and /„„ = C^ u . From the previous treatment it is clear that this amplitude again does not lead to any 
corrections of 0(T 2 ). 
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